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[t I ' The large-scale structure of the cosmos provides the key to understanding the early evolution of our Universe. 
I \ In particular primordial density perturbations may be produced from quantum fluctuations during an inflationary 
■ era. Most gauge-invariant studies of linear perturbations about homogeneous and isotropic cosmological background 
solution follow from the seminal article of Bardeen 0. In particular the reviews of Kodama and Sasaki ^ and 
psj . Mukhanov, Feldman and Brandenberger 3] have become widely used reference works. An alternative "covariant" 
J> ' approach has been advocated by Ellis and collaborators Q and developed by Challinor and Lasenby . These works 
, have provided gauge-invariant descriptions of the evolution of the total density, pressure and velocity perturbations 
in various cosmological models. 

Recently there has been increasing emphasis on the evolution of perturbations in multi-component systems where 
I it is useful to identify gauge-invariant adiabatic and entropy modes. This is partly driven by more detailed models 
of the multi-component systems in the early universe, such as reheating at the end of inflation 6], or the effect of 
. late-decaying scalar fields 7] . But it is also needed to make more detailed comparison with observational data in the 
late-time universe containing radiation, matter and dark energy (see, e.g., j^). 

On large scales, where the perturbed universe can be treated as locally homogeneous and isotropic, we can identify 
I a conserved perturbation with any locally conserved quantity [^. In particular there is a conserved perturbation 
O ■ associated with any perfect fluid whose energy density is locally conserved j^. This can be generalised to adiabatic 
perturbations of a multi-component cosmology and where the conserved perturbation is the gauge-invariant curvature 
perturbation on uniform total density hypersurfaces, introduced by Bardeen and collaborators |lllll2j |. 

Conversely, C, is not in general conserved on large scales when one considers non-adiabatic perturbations of a multi- 
component system. Nonetheless, at any instant, it is possible to decompose arbitrary perturbations of the multi- 
component cosmology into an adiabatic mode, with an associated conserved perturbation on large scales, and entropy 
^ ' (or isocurvature) modes which evolve independently of the adiabatic mode on large scales. One consequence of the 
5^ , effect of non-adiabatic perturbations upon the evolution of the total curvature perturbation, Q, can be the generation 
of correlations between the primordial curvature perturbation and any surviving isocurvature perturbations at the 
time of last-scattering of the cosmic microwave background . 

In this paper we will develop the gauge-invariant theory of cosmological perturbations to deal with multiple inter- 
acting fluids and fields. We develop the multi-component formalism of Kodama and Sasaki Q (see also 14] ). We 
subtly redefine the relative entropy perturbations to ensure gauge-invariance when including ener gy t ransfer between 
fluids. The resulting simplified evolution equations in the large-scale limit were presented in Ref. ^M,- We also relate 
this fluid description to previous work identifying adiabatic and entropy modes in multiple scalar field systems [T^ 
(see also Refs. jl^ lisl ITol 120^. 21]). Hwang has also investigated the evolution of perturbations with multiple fluids 
|22|. See Refs. [23 for an alternative long- wavelength approximation. 

The paper is organised as follows. In Section ^ we introduce the governing equations with multiple interacting 
fluids in the homogeneous background and for first-order inhomogeneous perturbations on all scales. In Section lTTTl we 
identify gauge-invariant perturbations describing the total density, velocity and shear perturbations and the relative 
entropy and velocity perturbations between different components, allowing for energy and momentum transfer between 
the fluids. We give the coupled first-order evolution equations and in particular give the simplified set of evolution 
equations valid in the large-scale limit. Including interactions enables us to treat fluids and scalar fields on the same 
footing and we show in Section Hvl how N scalar fields can be described as N "kinetic fluids" interacting with one 
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"potential fluid" describing the vacuum energy. We conclude in the final section with a short discussion. 

In Appendix^ we present the gauge transformation properties of the metric and matter variables including energy 
and momentum transfer and give some identities relating different variables that prove useful in simplifying our 
equations. In Appendix^ we give the governing equations on all scales in the uniform curvature gauge which may 
be better suited to numerical solutions. 

Throughout this work we assume a flat Friedmann-Robertson- Walker (FRW) background spacetime. Greek indices, 
/i, A, run from 0, ... 3, while lower case Latin indices, i, j, fc, run from 1, ... 3. Greek indices from the beginning of 
the alphabet, a,/?, 7 will be used to denote different fluids, and upper case Latin indices, I,J,K, denote different 
scalar fields. We only consider linear scalar perturbations in this article and leave vector and tensor perturbations to 
a future publication [23| . 



II. GOVERNING EQUATIONS 



In this section we give the governing equations for a system of multiple interacting fluids. 
The covariant Einstein equations are given by 

G fii, — SttG Tfj^iy , (2.1) 

where G^i, is the Einstein tensor, T^^ is the total energy-momentum tensor, and G is Newton's constant. Through 
the Bianchi identities, the field equations (|2.1(l imply the local conservation the total energy and momentum, 

(2.2) 

In the multiple fluid case the total energy-momentum tensor is the sum of the energy-momentum tensors of the 
individual fluids 

T^^-T.^t:)- (2.3) 

a 

For each fluid we write the local energy-momentum transfer 4-vector as 

"^.T^:) = QU > (2-4) 

where energy-momentum is locally conserved, Q'^^^) ~ ^- °^^^y ^'-'^ non-interacting fluids. Equations (|2.2I) and l|2.4(l 
imply the constraint 

E^M^O- (2.5) 

a 



A. Background equations 



The Einstein equations (|2.1|l give the Friedmann constraint and evolution equation for the background FRW universe 

= —P, (2.6) 

H = -4TrG{p + P) , (2.7) 
and energy-momentum conservation, Eq. H2.2|l . gives the continuity equation 

p=-'SH{p + P), (2.8) 

where p and P are the total energy density and the total pressure, a dot denotes a derivative with respect to coordinate 
time, t, the scale factor is a = a(t), and H = d/a is the Hubble parameter. 

The total density and the total pressure are related to the density and pressure of the component fluids by 
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The continuity equation (|2.4() for each individual fluid in the background is 

p„ = -3i/(p„ + P„) + Qa, (2.10) 

where the energy transfer to the a-fluid is given by the time component of the energy-momentum transfer vector 
Qfa) ~ ™ background. Equation H2.5|l implies that the background energy transfer obeys the constraint 

j2Qo. = o. (2.11) 



B. Perturbed equations 



We will consider linear scalar perturbations about a spatially-flat FRW background model, defined by the line 
element 

ds^ = -(1 + 2(t))dt^ + 2aB^,dtdx' + [(1 - 2V')<5y -I- 2E^,j\ dx^dx' , (2.12) 

where we use the notation of Ref. Q for the gauge-dependent curvature perturbation, ^, the lapse function, 0, and 
scalar shear, CTs = a?E — aB, and where 5ij denotes the flat background metric and X^i = dX/dx^ . The transformation 
of these scalar metric perturbations under arbitrary gauge-transformations is given in Appendix lA II 



1. Energy and momentum conservation 
The perturbed energy transfer 4- vector, Eq. 12. 4() . including terms up to first order, is written as llRj 

Q(a)0 = -(3a(l + 0) - <5Qq , 
Q{a)i = {fa + QaV)^i , 

and Eq. (|2.5|l implies that the perturbed energy and momentum transfer obey the constraints 



(2.13) 



(2.14) 



The perturbed energy conservation equation for a particular fluid, including energy transfer j_2S then obtained by the 
first-order part of the time-component of the perturbed continuity equation H2.4|l to give |2, |2^ 



Sp^ + 3H{Spo. + SP^) - 3 (/7„ + Pa) i! + —{po. + Pa) {Va + CTs) - Qa4> - <5Q„ = , 



(2.15) 



where 6pa and 5Pa are the perturbed energy density and the perturbed pressure of the a-fluid, respectively, and the 
comoving spatial Laplacian is denoted by = S"^^ /dx^dxL Va is the covariant velocity perturbation of the a-fluid 
defined as 

Va = a{va, + B), (2.16) 
where Va is the scalar velocity potential of the a-fluid. The momentum conservation equation of the a-fluid is 



Qa 



Pa+Pc 



il + cl)~3Hcl 



1 



Pa + Pa 



2 V2 

SPa + ^^Ila-QaV-fa 



= 0, 



(2.17) 



where IIq, is the perturbed scalar anisotropic stress of the a-fluid and = Pa/ Pa is the adiabatic sound speed of that 
fluid. Note that for stress-free, non-interacting dust we have Va + 4> ~ and hence Va =constant in a synchronous 
gauge, whereas the scalar velocity potential redshifts as Uq oc 1/a. In Ref. j3| we used 5qa = {pa + Pa)Va as the 
momentum perturbation. 

The total fluid perturbations are related to the individual fluid quantities by 



Y,6pa^6p, Y.^Pa = 5P, ^n„ = n, 

<x a. OL 



(2.18) 
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and 
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where 5p is the total density perturbation, 5P the total pressure perturbation, 11 the total anisotropic stress, and the 
total covariant velocity perturbation is given by 

V = a{v + B), (2.20) 

where v is the total scalar velocity potential. 

We therefore get an evolution equation for the total density perturbation from Eq. H2.15() by summing over all 
fluids, using Eq. H2.18|l and the constraint H2.14|l . 

• V2 

5p + ?,H{5p + ,5P) - 3 (p + P) + ^(p + P) (1/ + (Is) = , (2.21) 

while the total momentum conservation equation is given from Eq. H2.17|l . using Eqs. l|TT5|l and l|TT^ and the 
constraint (|2.14|l . 

V -iHclV + cl,+ ^-(5P+\%n]=Q, (2.22) 
p + P \ 3 / 

where Cg is the adiabatic speed of sound, defined as 



< = -■ (2.23) 



2. Einstein's field equations 

The Gq component of the Einstein equations H2.1|l yields the first-order perturbed energy constraint equation 0, 0, 

3iJ + - + -f^f^s) + 47rG(5p = . (2.24) 

We get the momentum constraint equation (identically zero in the FRW background) from the G°-component of the 
Einstein equations, 

+ H</> + 47rG(p + P)y = . (2.25) 

From the trace-free part of the component of Einstein's equations we find the shear evolution equation 

6^ + Ha^ - (/) + -0 - SttGH = , (2.26) 

and the trace of the G* component of Einstein's equations gives, using Eq. ()2.26() . the evolution equation for the 
curvature perturbation as 

^ + + H4:+ {iH^ + 2H^ 4> - 47rG (sP +'^^^ = . (2.27) 



III. GAUGE-INVARIANT PERTURBATIONS AND GOVERNING EQUATIONS 

In this section we define gauge-invariant variables, using the transformation properties of the metric and matter 
perturbations given in Appendix lA II We then use the results of Section III BI to get the evolution equations for these 
gauge-invariant variables. 

In this and the following sections we will replace spatial Laplacians by the wavenumbers of their respective eigen- 
modes according to — > —k^. 
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A. Definitions 



1. Curvature perturbations 



A gauge-invariant definition of the curvature perturbation for each fluid is given by [l£ 

Ca = -1P - H^- . 

Pa 



(3.1) 



This describes the dimensionless density perturbation on uniform-curvature hypersurfaces, or, equivalently at first- 
order, the curvature perturbation on uniform a- fluid density hypersurfaces |2fi| |. The total curvature perturbation on 
uniform density hypersurfaces [lol Hit IT^ 



(: = -4,-H 



Sp 



(3.2) 



is then a weighted sum of the individual fluid perturbations 

p 



(3.3) 



We can also deflnc a gauge-invariant combination of the curvature and velocity perturbations for each fluid 

na=Tp'- HVa ■ (3.4) 

This describes the curvature perturbation on hypersurfaces orthogonal to worldlines comoving with the fluid, and is 
particularly useful when studying scalar field perturbations. The total comoving curvature perturbation is given by 

MM 



and is again a weighted sum, using Eq. H2.19|) , of the individual fiuid perturbations 



Pa 



p+p 



Ha- 



(3.5) 



(3.6) 



The two alternative descriptions of the total curvature perturbation are closely related. Indeed even the coefficients 
Pa/p appearing in Eq. H3.3|l and (pa + Pa) / {p + P) in Eq. H3.6|l are the same for non-interacting fluids, i.e., when 
Qa = 0. The constraint equation H2.24|l can be rewritten in terms of these gauge-invariant quantities as 



(3.7) 



where we introduce the gauge- invariant curvature perturbation in the longitudinal or zero-shear gauge P, Q , 

* = V + iJcTs . (3.8) 

In practise it is convenient to use only one of C or 7?. and use 5* to keep track of C + which, as can be seen from 
(|3.7|) . typically becomes small on super-Hubble scales, i.e., k^/a^H^ ^ 1. 
The evolution equation for (a is given from Eq. (|2.15|) and using Eq. (|2.24|) 



Ca 



3— iSPa " C^bpa ) - H — 
Pa ' Pc 



H 

Pa 



SQa - ^Spo 
Pa 



* - ( 1 - ^ ) 7^a 

Pa 



(3.9) 



where the terms in the first line are gauge-invariant combinations, which we shall deal with in the next subsection on 
entropy perturbations. The evolution equation for TZa is given from the momentum conservation equation, Eq. 12.17|l . 
and using Eq. 



Tla = ^(7^„-7^) ^ 

H pa + 



Pa 



Hcl (7^„ + Ca 



H 



Pa + Pa 



5Pa - clSp^ - ^^n, -fa-QaiV~ V^) 

3 a-^ 



(3.10) 
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Finally we note that the curvature perturbations on uniform total density, uniform shear and comoving hypersur- 
faces, ^, ^ and TZ, are related to the density, shear and the velocity perturbations on uniform curvature hypersurfaces, 

and similarly for the curvature perturbations on uniform a-density and uniform a-velocity hypersurfaces 

where the "tilde" denotes perturbations on uniform curvature hypersurfaces. We can therefore also think of the cur- 
vature perturbations C and (a as describing density perturbations or density contrasts, of the curvature perturbations 
TZ and TZa as describing velocities, and of the curvature perturbation ^/ as describing the shear. 

Governing equations for the variables defined in Eqs. (|3.11() and and (|3.12l) in uniform curvature gauge are given 
in Appendix IBl 

2. Entropy perturbations 

The difference between the density perturbation for any two fluids describes a gauge-invariant relative entropy 
(or isocurvature) perturbation 0, ^3 

Sc.f3 ^ 3(Ca - C/3) = -3H _ (3.13) 

\ Pa PfJ J 

Note that this only coincides with the entropy perturbation defined by Kodama and Sasaki Q in the absence of 
energy transfer, Qa = Qp = 0. The factor of 3 is introduced to coincide with the conventional definition of entropy 
perturbations in the baryon-photon ratio: 

'-^^ = 5b, , (3.14) 

nB/n, 

where jt-b and n-y are the baryon and the photon number densities, respectively. 
Similarly we define the gauge-invariant relative velocity perturbation by |2| ^ 

V^p = V^-Vp^-^ (7^„ - TZp) . (3.16) 

The total pressure perturbation can be split into an adiabatic and non-adiabatic part 

5P = 6P^^d + cl6p, (3.17) 

where Cg is the adiabatic speed of sound defined in Eq. H2.23|l above. 

In the presence of more than one fluid, the total non-adiabatic pressure perturbation, ^Pnad, may be further split 
into two parts 0, 

5P„ad = <JPi„tr + <JProl . (3.18) 

The first part is due to the intrinsic entropy perturbation of each fluid 

<5Pi„tr = 51 '^^intr.a , (3.19) 



^ The definition of the relative velocity perturbation in is related to our definition in Eq. 13.161 by 



(3.15) 
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where the intrinsic non-adiabatic pressure perturbation of each fluid is given by 



(3.20) 



For any fluid with a deflnite equation of state, Pa = Pa{Pa), the intrinsic non-adiabatic pressure perturbation must 
vanish, ^Pjntr.a = 0. The total adiabatic sound speed 12.23|l is the weighted sum of the adiabatic sound speeds of the 
individual fluids, 



^ P 



Pa 2 



(3.21) 



The second part of the non-adiabatic pressure perturbation H3.18|l is due to the relative entropy perturbation Sap 
between different fluids, deflned in Eq. H3.13|l . 



a,/3 



^a. 



P ■ 



(3.22) 



Analogous to the non-adiabatic pressure perturbation for each fluid H3.20|l . we can identify an intrinsic non-adiabatic 
part of the energy transfer perturbation 15] that appears in the perturbed continuity equation for each fluid (|2.15(1 



5Qa 



-Spa ■ 



(3.23) 



This is automatically zero if the local energy transfer Qa is a function of the local density pa so that SQa = 
(Qa/ Pa)Spai just as the intrinsic non-adiabatic pressure perturbation (|3.2U|) vanishes when SPa = {Pa/ Pa)^Pa- Note 
however, that from the deflnition, Eq. H3.23|l . the intrinsic non-adiabatic energy transfer perturbation, ^Qmtr.a, can 
be non-zero in the case where the background energy transfer Qa ^ 0, even if 5Qa = 0. We can also identify a relative 
non-adiabatic energy transfer that appears in Eq. (|2.15|l whenever Qa 7^ 



AH — f ^ 

H \ Pa P 



Qa 



•af3 : 



(3.24) 



due to the presence of relative entropy perturbations. Only the intrinsic parts of the non-adiabatic pressure and 
energy-transfer perturbations will appear in the evolution equation for the relative entropy perturbation, Eq. H3.31|) 
below. 

With the above definitions the evolution equation for the curvature perturbation Cq, Eq. H3.9|) . reduces to 



H 



1 fc2 



3-; — 5Pintr, Q ;— {5Q\ntr,a + 5Qrc\,a) + 7777^ 

Pa Pa in a- 



* - ( 1 - — ) 

Pa 



(3.25) 



Thus we see that ~ in the absence of non-adiabatic pressure and energy transfer perturbations, on the large 
scales where we can neglect gradient terms. 

There is no momentum transfer in the background FRW universe so the momentum transfer fa that appears in the 
momentum conservation equation for each fluid (|2.17|) is automatically gauge-invariant. Nevertheless we can identify 
a relative momentum transfer in Eq. (|2.17l) . defined as 



/rcl,a = Qa {V — Va) , 

which can be rewritten in terms of the relative velocity perturbation, defined in Eq. H3.16|l . using Eq. 



(3.26) 



2l, as 



Ac 



EPj_+J^ 



p + p 



The evolution equation for TZa , Eq. H3.10(l then simplifies to 



^a = J^ {na - TZ) 



Pa 



Pa 



Hcl {TZa + Ca) 



H 

Pa + Pa 



2 

^Pintr.a 77 0-^0 fa /rcl.a 
3 



(3.27) 



(3.28) 



We end up with the following set of gauge-invariant dynamical variables: the curvature perturbation on uniform den- 
sity hypersurfaces, C, the curvature perturbation on uniform shear hypersurfaces ^f, the relative entropy perturbation 
SafSi and the relative velocity perturbation Vap- 
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B. Evolution equations 



In this sub-section we re-express the gauge dependent evolution equations given in Section III Bl in terms of the 
gauge- invariant quantities defined above in Section [ill Al In deriving the evolution equations, we make extensive use 
of the relation between the relative entropy perturbation and the a-fluid and total density perturbations, Eq. |IA23|) . 
and the relation between the relative velocity perturbation and the a-fluid and total velocity perturbations, Eq. ljA24p . 

The first-order evolution equation for the total curvature perturbation, ^, defined in Eq. 13. 2|) . derives from the 
evolution equation of the total density perturbation, Eq. H2.21(l . and is given by 



c 



H 



P 



nad 



1 fc2 
fljf 



C) 



1 



9HH 



(3.29) 

defined in Eq. 



where SPnad, as defined in Eq. (|3.18|l depends on the intrinsic entropy perturbations, SPir 
and the relative entropy perturbations, Sa/s defined in Eq. H3.13|l . 

The evolution equation for the curvature perturbation in the zero-shear gauge, ^P, defined in Eq. H3.8|l . follows from 
Eqs. (j^^ and fT^H)! 



H - 



H 
H 



1 P 
Iff ^ 



H 
H 



C = -SttGHU . 



(3.30) 



Note that energy-momentum transfer between individual fluids does not directly enter into the evolution equations 
for the total curvature perturbations C and 5*. However it does considerably complicate the evolution equations for 
the relative entropy and relative velocity perturbations, as can be seen below. 

The evolution of the relative entropy perturbation. Sap defined in Eq. H3.13|l . using Eq. (|2.15|) . is given by 



Saf3 = 



fSHdPi 



intr,Q 



/intr 



3HSP 



intr,/3 



/inti\/3 



H 
2H 



-. H Oq/3 




(K.7 + Va^) 



(3.31) 



Only the intrinsic non-adiabatic perturbations 5-Pintr.Q and 5(5intr,Q appear explicitly in the evolution equation H3.31|l . 
while the relative part of the non-adiabatic energy transfer perturbations H3.24|l have been expressed in terms of Sap. 

Finally, using the evolution equation for the velocity perturbation of the a-fiuid, Eq. (|2.17|) . we can derive an 
evolution equation for the gauge-invariant relative velocity perturbation, Vap defined in Eq. 13.16|l . 



a/3 




(3.32) 

-TTa fjp-tJrpy 

where we have used Eq. (j3.27|l for the relative momentum transfer perturbation. 

We therefore see that the evolution equations for (, Sap and Vap, Eqs. (|3.29l) - (|3.32|) . form a system of coupled, 
first order, ordinary differential equations. Although we have defined n{n — l)/2 relative entropy perturbations, Sap, 
for n fluids, there are of course only n — 1 independent entropy perturbations, and similarly n — 1 independent relative 
velocity perturbations. Thus we have 2n coupled equations relating 2n variables describing the density and velocity 
perturbations in n fluids. 
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The system is not closed until we specify the intrinsic non-adiabatic pressure perturbation 5Pintr,a, the intrinsic non- 
adiabatic energy transfer perturbation (5(5intr,a) the anisotropic stresses IIq, and the perturbed momentum transfer /q, 
for each fluid. In the background we have to prescribe the equation of state, Pa{pa), which determines the adiabatic 
sound speed c^, and the energy transfer Qa- All these quantities depend upon the physical model for the fluids. For 
perfect fluids 5Pintr,a, SQintr,a, Hq, , and fa are all zero. In Section Hvl we will calculate these quantities explicitly for 
the particular example of N scalar fields. 



C. Large-scale limit 

Throughout much of the early history of our Universe, scales of astrophysical interest are far larger than the Hubble 
scale, which defines the cosmological expansion time. For instance, the Hubble scale at the epoch of primordial 

nucleosynthesis corresponds to a comoving scale of around 10 pc in the present Universe. Thus it is often a very good 
approximation to work in a large scale limit when studying the origin and evolution in the early universe of the large 
scale structure of our observable Universe. In this section we give the governing equations in the large scale limit, 
i.e. in the limit fc — > 0. Some of the results in this section have already been presented in |15|. 

The form of the evolution equations for the perturbations suggests that the time-dependent amplitude of any 
perturbation, 5X{t)e^^'^ , can be expanded as a Taylor series in terms of the comoving wavenumber, fc, as 



5X{t) 



ri,=0,l,... 



(3.33) 



Equations (|3.29l) - (|3.32l) then give coupled first-order evolution equations for X„(i) driven by X„_i(t) for n > 1 [and 
driven by Xn-2 for n > 2 in Eqs. and ^^]. 



The long- wavelength limit as k 
of first-order equations: 



is given by the time-dependence of the n — Q solutions to the homogeneous set 



o(0) 



V, 



afi 



H 



H 



(0) 



= 3iJ 



(0) _ 



Qo 



Pa + Pa 

Q a 
Pa + Pa 



Pf> ■ 



H 
H 



-,(0) 

■' intr,a 



6Q 



(0) ^ 
intr,/3 



Pfi 



(0) 



a/3 



(1 



P/3 



1 

' 2 
-(1 



P/3 



E 



P7 f Qa 



2p \ Pa 



Pa + Pa PI3 + Pf} 



3H {cl 



Pa 

„2 ^ 



Pa PP 



Pr 



15(0) 



pp 



r(0) 



(3.34) 
(3.35) 

(3.36) 



'iH id + cl\ 



Vap - —O, 



(0) 

^JJ^aP 



•yy Pj_ 

7 



P. 



P 



fy(0) 
\^ ai 



V, 



(0) 



Pi 



+ 



AO) 

Ja 



f(0) 
Jp 



sp: 



(0) 



5P 



7 

(0) 



Pa PP+ Pp 



Pa PP + Pp 



(3.37) 



This long-wavelength limit exists so long as the inhomogeneous source terms in the full evolution equations 
(|3.32l) . that is i5Pintr,Q, ^Qintr.a, fa and Ha, all have a well-defined limit as A: — > 0. 

The perturbation, ^, vanishes from the other evolution equations completely in the large scale limit. Since this 
variable describes the shear (on uniform curvature hypersurfaces), its disappearance can be intuitively understood by 
observing that the shear is expected to vanish on large scales. Because the shear and velocity perturbations, ^' and Vap, 
only enter the evolution equations (|3.29|l and H3.31|l at order k^, we can obtain a long- wavelength limit for for C^"-' and 

intr,Q, k'^fa and k'^Ua all have a well-defined hmit as A: — > 0. This 



^ap under the weaker requirement that 6Pintr,a, 



is the basis of the 'separate universes' picture jlOl | commonly used to study the evolution of density perturbations on 
sufficiently large scales where the universe looks locally like an unperturbed (FRW) cosmology. Specifically it requires 
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that we can neglect the divergence of the momenta in the zero-shear gauge, V^[(/3q + Pa){Va + cts)], in the perturbed 
continuity equation H2.15|l for each fluid. 

We see from Eq. H3.34|l that in the large scale limit the evolution for the curvature perturbation C, is only sourced 
by the non-adiabatic pressure perturbation, as defined in Eq. (|3.18|) . including contributions from the intrinsic non- 
adiabatic pressure perturbations of each fluid and the relative entropy perturbations between fluids. For dPn^d — 
we recover the famous result that ( =constant for adiabatic perturbations on large scales |3, Hfl Ull Il2 . 12M l29l | . The 
long-wavelength solution for adiabatic density perturbations is 



(3.38) 

where C and D are constants of integration. This adiabatic solution exists on large scales even in the presence of 
energy transfer between fluids, Qa 7^ 0, so long as the intrinsic non-adiabatic energy transfer, ^Qintr.a defined in 
Eq. H3.23|l . is zero for each individual fiuid. 

The evolution of the relative entropy perturbations, Sap, is independent of the curvature perturbations C and 
'3/ in the large-scale limit, Eq. H3.36I) . The actual evolution is dependent upon the intrinsic non-adiabatic pressure 
and energy transfer of each fluid. However for perfect fluids all the source terms vanish we have the simple result 
Sa/3 =constant in the large-scale limit. 

The evolution of the relative velocity perturbations Vaf3, Eq. H3.37|l is driven by relative entropy perturbations, Saf3, 
as well as the intrinsic non-adiabatic pressure perturbations and momentum transfer, even on large scales. Only for 
pressureless, non-interacting fluids do we find Vaf3 =constant on large scales. 









H 




a 


c(0) 


= 0. 



D 



~ 8^Gn(") j dt 



IV. SCALAR FIELDS 



In this section we show how scalar fields can be included in our interacting fiuid formalism. Multiple scalar fields 
provide a test case in which quantities like the energy transfer and the relative momentum transfer are simple enough 
to be computed from first principles, and can then be applied to identify adiabatic and entropy modes in the fluid 
formalism. It also enables us to later "mix" fluids and scalar flelds in a consistent way, which is necessary in many 
cosmological models, including quintessence models of dark energy, decaying scalar fields, such as the curvaton, or 
inhomogeneous reheating at the end of inflation. In particular we will identify the adiabatic mode and relative entropy 
perturbations in the fluid description and compare these with previous analyses of adiabatic and entropy modes in 
the scalar held perturbations [ig . 

We will split the total energy-momentum of N scalar fields into N maximally stiff fiuids interacting with a potential 
energy. Our splitting is slightly different from that of Hwang and Noh who also relate the scalar field perturbations 
to fluid quantities. They work in terms of the total density and pressure perturbation for each field, including kinetic 
and potential perturbations for each field. There is thus an intrinsic non-adiabatic perturbation for each field. By 
contrast we separate the kinetic and potential perturbations so that each has a fixed equation of state and thus no 
intrinsic non-adiabatic pressure perturbation and we can always relate non-adiabatic perturbations to relative entropy 
perturbations between different component parts. 

For simplicity we consider here only minimally coupled scalar fields in the sense that they interact only via their 
combined self-interaction potential. Although we allow for the presence of fluids other than the scalar fields, we 
assume that the scalar fields do not exchange energy-momentum with any other fiuids. 
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A. Background 

1. Standard treatment 

The energy density and the pressure of N minimally coupled scalar fields, labelled by the subscript are at 
zeroth order 



= ^E^?-^' (4.2) 



2 



where U = U((pi) is the potential of the scalar fields and upper case Latin indices, /, J, K, run from 1, . . . , iV. 
The adiabatic sound speed of all the scalar fields is given from Eq. 12.23|l by 



where U = J2i U^^jipi. 



2. Splitting into kinetic and potential parts 



We can split the total energy density and pressure of N scalar fields into N "kinetic fluids" with energy density 
and pressure, respectively, 

PI = \v] , Pi = , (4.4) 
and a single "potential fluid" or vacuum energy, 

Pu = U, Pu = ~U. (4.5) 

We now have 

The kinetic and potential fluids have a barotropic equation of state and adiabatic sound speeds 

Cs/ = 1, clu = -l. (4.7) 

The Klein-Gordon equation for each field, 

(pi + iHi^i + U,^, ^Q, (4.8) 

gives a continuity equation for each kinetic energy H4.5() of the form given in Eq. where the energy transfer to 

the kinetic fiuids from the potential is given by 

Qi = --ipiU,^, . (4.9) 
where U^^j = dU/dLpj. The constraint ()2.5|l then implies that the energy transfer to the potential is 

Qc; = ^(^/C/,^,. (4.10) 
I 

B. Perturbations 



We now extend the fluid formalism to the perturbed scalar field case. 
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1. Standard treatment 



The total energy density and the pressure perturbation for N scalar fields are 



(4.11) 
(4.12) 



having contributions both from the individual fields and the potential, whereas the total velocity perturbation, given 

by 



(4.13) 



is independent of the potential. Scalar fields cannot support anisotropic stress, and we have 

n^ = 0. (4.14) 
The total non-adiabatic pressure perturbation for N scalar fields is then readily calculated from Eq. (|3.17|l to be 

2 



6Pr 



</3/ ^ipiSU - U {dip J ~ ipicjjj 



where 



(4.15) 



(4.16) 



by 



The total entropy perturbation for N scalar fields was denoted S in Ref. [iQI ^ind is related to the above definition 



H 



The total entropy perturbation reduces in the single field case to 

2C/,^(p 



SR 



nadc^ 



3H 



dip — ip<j) Sip 
p ip 



which is related to the intrinsic entropy perturbation T defined in Ref. |30i by 
and to the intrinsic entropy perturbation Fq defined in Ref. |3l| by SPn^d^p = ^VTq. 



(4.17) 



(4.18) 



(4.19) 



2. Splitting into kinetic and potential parts 

Splitting the energy and pressure perturbation, as for the background fields, into N "kinetic fluid" quantities and 
a single "potential fluid" we get 



5pi = ipiSpi - ip](j), 
SPi = ipiSpj - ipj(f), 
Hi = 0, 



(4.20) 
(4.21) 
(4.22) 



for each kinetic fluid, and 



Spu = SU ^'^U^^jSpj , SPu = -'^U^^jSpj , Uu^Q, 



(4.23) 



13 



for the potential energy. We thus have N gauge-invariant curvature perturbations, C/ defined according to Eq. (|3.1|l . 
for each kinetic fluid plus one gauge-invariant curvature perturbation, for the potential energy. 
The covariant velocity perturbation for each kinetic fluid is 

Fz = -^. (4.24) 

Substituting Eq. (|4.24|) into the definition of the total velocity, Eq. 1)2. 19|) . we get the standard result, Eq. (|4.13l) . for 
the total velocity perturbation due to the scalar fields. Note that the potential energy has no momentum perturbation 
and hence its covariant velocity perturbation, Vu, is undefined. 

The comoving curvature perturbation, defined in Eq. (|3.4|) . for each field is 

TZi = iP + H^, (4.25) 

and hence the total comoving curvature perturbation (i.e. relative to the average fluid velocity) is given by 

The Sasaki-Mukhanov variable j33. l33l| . or the field fluctuation on miiform curvature hypcrsurfaces, is simply Qj = 

Note that the held perturbations determine not only the velocity perturbations (|4.24(l but also the potential energy 
perturbation H4.23|l . Thus 5pu and V^p are not independent variables. In terms of the gauge- invariant comoving 
curvature perturbations TZi we then have, using Eq. (|A26(I . 



Cu ^-j Y^^^vk^kTIk, (4.27) 

which we can rewrite in terms of gauge-invariant total scalar field comoving curvature perturbation H4.26|l and velocity 
perturbations (|4.13|l and H4.24|l 

Qu = -n^ + HY %^ {Vj -V^). (4.28) 

For a single scalar field we have simply C^u = —TZ^. 

Although we have decomposed the total energy-momentum of the N scalar fields into N kinetic fiuids interacting 
with one potential "fluid" , the potential has no momentum perturbation and its density perturbation can be written 
in terms of the kinetic fluid velocities. Hence we only need the TV density perturbations and N velocity perturbations 
to describe the perturbed energy-momentum of N scalar flelds. 

We get the perturbed energy transfer to each field by substituting Eq. (|4.20|l for the fluid density into the perturbed 
energy conservation equation 1)2. 15|) . and using the perturbed Klein-Gordon equation for each field, 

U2 



5(fj + 3H5(pj + — {5lpj - (fiias) + ^ U^^,^j5(fj -ipi(^^ + Sipj + 2J7,^,0 = , (4.29) 
" J 

which gives 

SQj = -f/,<pj (Sipj - (y5/(/>j - 95/ X! ^.vivj^^j ■ (4.30) 

J 

The perturbed energy transfer constraint, H2.14|l . then gives the energy transfer perturbation to the potential 

6Qu =5U -U(j), (4.31) 

where 

5U = Y [^iU,^,^M.] + E ^^^'^'^i ■ (4.32) 

I 
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Finally we identify the momentum transfer to each component. Substituting the expressions for the /-fluid quan- 
tities, Vj, 6Pj, pj + Pj and cj of Section llV Ij 21 into the momentum conservation equation H2.17II and, noting that 
Hi — 0, we obtain the momentum transfer perturbation to the /-fluid 

// = t/.^x {S^i + ^iV) , (4.33) 

in terms of the field variables or, using Eq. H3.26|l . 

// - Qi {Vi ~V)^ -Aci,i . (4.34) 

Hence the total momentum transfer for each field, // -I- /roi.i, that appears, for instance, in Eq. (|3.28|) . is zero. 

As noted earlier the potential fluid has vanishing momentum, and from the momentum conservation equation (|2.17|) 
for the potential we require 

fu = QuiV-V^), (4.35) 

which is zero if there are no fluids present other than the scalar fields and V = V^p. One can verify that Eq. 12.14(1 for 
the total momentum transfer is satisfied with fu + J2i fi = 0- 



C. Relative perturbations 



Substituting the expression (|4.2UI) for the density perturbation of each "kinetic fluid" , and the background density 
(14.4(1 . into the definition for the relative density perturbation, Eq. 1(3.13(1 . we get 



-3H 



If I </3,/ 



(4.36) 



which is, despite appearances, gauge-invariant as can be seen from the gauge transformation properties of the variables, 
Eqs. and ljIT8)l . 

The relative density perturbation between each field and the potential is given by 



Sep J — (picf) SU 
(pi u 



(4.37) 



However, remembering that the density perturbation for the potential can be written in terms of the field velocity 
perturbations, as in Eq. ((4.28() . we can also write the relative entropy perturbation, Eq. ((3.13(1 . between any fluid and 
the potential as 



= 3(Ca+7ev)-3/r^ 



u 



(4.38) 



In the presence of additional fluids the total comoving curvature perturbation 1(3.6(1 is given by 

Pa + Pcx 



p + p 



and if we then use Eq. 1(3.7(1 to eliminate 7?. in favour of C and we finally obtain 



Sal 



p-U 



1 fc^ T P7 n sr^ TT . Pi + P-1 



(4.39) 



(4.40) 



where we also used Eq. ((A25(l . 

The definition ((4.36(1 of the relative entropy perturbation, 5/7, between two fields in terms of their relative density 
perturbations, 5pi, differs from the definition used in Refs. [liL l20l l2ll| . where the relative entropy perturbation 
was defined in terms of the relative field perturbations, 5lpj. In our fluid formalism the relative field perturbation is 
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proportional to the relative velocity perturbation, defined in Eq. (|3.16|l . which in the scalar field case, using Eq. H4.24(l . 
is given by 

Vu = -('-^-'-^). (4.41) 

Indeed we can express the relative velocity perturbation, defined in Eq. 1)4.41(1 . as the difference between comoving 
curvature perturbations, 

Vi., = -^[ni-nj) , (4.42) 

which is analogous to the definition of the relative entropy perturbation Sap expressed as the difference of curvature 
perturbations on uniform density hypersurfaces in Eq. I|3.13|l . We shall see in a moment that when we evaluate the 
non-adiabatic energy transfer in terms of the relative perturbations between fields, the relative velocity perturbation 
(|4.41() will indeed appear as a non-adiabatic perturbation. 

Note that, as Vjj is not defined for the potential energy, the relative velocity perturbation, Vjja, between any 
fluid and the potential energy is not defined. This would be a problem if Vjja appeared in the evolution equations 
(|3. 291 153^ for our other fluid perturbations. However, Vua only appears in summations over all fluids, 7, where it is 
multiplied by pj + P^, which vanishes for the potential energy. 



D. Non-adiabatic perturbations 



Each kinetic fluid and the potential fluid have definite equations of state {6 Pi = dpi, dPjj 
(|3.2U|) . there is no intrinsic non-adiabatic pressure perturbation 



SP 



intr.I 



0. 



6Pi 



intr.U 



0, 



-Spu)- Hence, from 



(4.43) 



And since the kinetic fluids all have the same adiabatic sound speed {c^j = 1) there is no relative non-adiabatic 
pressure perturbation (|3.22|) due to entropy perturbations between the different kinetic fluids, Sij. The total entropy 
perturbation, Eq. (|4.15l) . is thus due solely to the relative entropy perturbation between the fields and the potential 
(which has c^^ = — 1): 



SP 



2U 



— X! V'lV'iSiu ■ 



(4.44) 



Note, however, that this entropy perturbation can always be re-expressed in terms of a sum over the other entropy 
perturbations and relative velocity perturbations, Eq. H4.38|l . 

The intrinsic non-adiabatic energy transfer perturbation, defined in Eq. H3.23|l . for the /-fluids follows from 
Eqs. iHini), lEHni), and 10201, and is given by 



^QintrJ = 'Pl 



P- (Sipj - ipict)) - SU,^j 



which can be rewritten as 



-Sip J + ^[Sipj 
ipi 



(4.45) 



(4.46) 



This is a relative entropy perturbation between the proper time derivative of the I-th field and the value of the J-th 
field. It is convenient to expand this as 



trJ = X! V'lV'jU^^j^j 



W_5^ {f^fii_^i^ _SU_ 

if ipj ipi u 



because, using Eq. 14.41(1 . we can write 

SU Sipj 



U 



SfK _ ^ 



IPK 



E 

K 



u 



VjK. 



(4.47) 



(4.48) 
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So finally, using Eq. (|4.37|) . we can express the non-adiabatic intrinsic energy transfer in terms of the relative density 
perturbations Sju and velocity perturbations Vjk 



. K 



(4.49) 



For the U-fluid the intrinsic non-adiabatic energy transfer perturbation is given by 



(^Ointr.u = 61J -U(l)- ^SU , (4.50) 
where we have used Eqs. H3.23|l . (|4.5I) . 1)4. 9|) . and (|4.23|) . which yields in terms of relative entropy perturbations 

SQintr,V = -T^ U^iViSiu + ^ U^,^,^I^ j^^^-^^^^^Vk J . (4.51) 

I I,J 

Thus we see that a non-zero relative velocity perturbation, Vjk , between scalar fields leads to a non-adiabatic energy 
transfer in Eqs. 14.49|l and 14.51|l . To recover the adiabatic solution 1)3.38(1 on long- wavelengths with Siu = we see, 
from Eq. (|3.36(l . that we require in addition that Vjk — *■ for scalar fields with arbitrary interaction potential U{ipi). 
Thus we can identify the relative velocity perturbation Vjk between scalar fields as a non-adiabatic perturbation 
^"6], in contrast to the case of perfect fluids where we require only that k'^Va/s ^ as A: ^ to obtain an adiabatic 
solution with Safj = on large scales. 

We can finally rewrite the momentum transfer H4.34|l solely in terms of a sum over the relative velocity perturbations 
Via, according to Eq. 13.27|l . 

// = -^/t/,,,^^|±^l/,„. (4.52) 



V. DISCUSSION AND CONCLUSIONS 

In this paper we have extended the analysis in the classic paper by Kodama and Sasaki to deal with the coupled 
evolution of curvature and isocurvature perturbations in a multi-component cosmology with interacting fluids and 
scalar fields. In doing so we have clarified the nature of adiabatic initial conditions, and identified non-adiabatic 
effects. 

It is well-known that the total curvature perturbation ^ is conserved on large scales when the total non-adiabatic 
pressure perturbation (|3.17|) vanishes. In a multi-component system the total non-adiabatic pressure perturbation 
can be split, according to Eq. H3.18|l . into the sum of the intrinsic non-adiabatic pressure perturbations of individual 
components and the relative pressure perturbation due to the relative density perturbations. Sap, between different 
components. The intrinsic non-adiabatic pressure perturbations depend on the internal degrees of freedom of the 
fluid, but must vanish for fluids with a definite equation of state PaiPa)- 

The relative density perturbation between different components, Sap defined in H3.13|l . generalises the perturbation 
in the photon-baryon ratio in a standard hot big bang cosmology to arbitrary interacting fluids. Hence we refer to 
Sap as a relative entropy (or isocurvature) perturbation. 

We have given, for the first time, the evolution equations on all scales for the gauge-invariant entropy perturba- 
tions Saf3 in a multi-component system, which is driven by the intrinsic non-adiabatic pressure perturbations and 
perturbed energy transfer. Analogously to the pressure perturbations we split the perturbed energy transfer into the 
gauge-invariant intrinsic non-adiabatic energy transfer and the relative energy transfer, due to the relative entropy 
perturbations, Sap- In certain models, such as spontaneous decay of non-relativistic particles [l5l |. one may also be 
able to give the intrinsic non-adiabatic energy transfer in terms of the relative entropy perturbations, enabling one to 
obtain a closed set of evolution equations. 

The overall curvature perturbation and relative entropy perturbations are also coupled to the velocity perturbations, 
but these decouple in the large scale limit (in the absence of intrinsic non-adiabatic pressure and intrinsic non-adiabatic 
energy transfer). Hence adiabatic density perturbations are characterised by a single amplitude, C, which remains 
constant on large scales. A sufficient condition for adiabatic perturbations on large scales is 

SPinti-,a = , (5Qi„tr,a = , and Sap = , (5.1) 
for all components a and (3. Adiabatic perturbations then stay adiabatic on large scales. 
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We have also shown how to describe scalar fields and their perturbations in this notation. N interacting scalar 
fields can be described by N kinetic "fluids" with a stiff equation of state, P/ = pj, interacting with one potential 
"fluid" with vacuum equation of state, Pjj = —pu- The potential fluid velocity is undefined, because its momentum 
perturbation {pjj + Pu)Vu is zero. Because the components have fixed equations of state there are no intrinsic non- 
adiabatic pressure perturbations, only the relative non-adiabatic pressure due to relative entropy perturbations. On 
the other hand there is an intrinsic non-adiabatic energy transfer which includes terms due to the relative velocity 
perturbation, Vjj, even on large scales. Thus adiabatic perturbations for interacting scalar fields on large scales 
require 

5/j = 0, 5/c/ = 0, and y/j = 0. (5.2) 

Each scalar field perturbation S^pi determines the velocity perturbation, Vi, and thus adiabatic field perturbations 
[Tfij l must obey 

Vij = — = . (5.3) 

During inflation scalar field perturbations originate as small-scale quantum fluctuations and these will not in general 
respect this adiabatic condition. On the other hand the existence of a unique attractor in phase space at late times 
would drive perturbations towards Vij = 0. In addition the adiabatic condition, Sij — 0, for relative density 
perturbations requires that 

^ = ^ , (5.4) 

VI VJ 

where 5lpit = Sipj — ipj(j) is the perturbed proper time derivative of the fleld, and Sju — requires 

Siplr SU 



VI u 



(5.5) 



All of the above conditions can be enforced by requiring the generalised adiabatic condition for perturbations in 
any 4-scalars x and y jlfl| : 

- ^ ^ . (5.6) 

X y 

This is a generalisation of the usual adiabatic condition, Saf} — 0, for fluid density perturbations which requires 

i_Pa_ ^^_P!i_ 
Pa pp ' 

Because the scalar field perturbation determines both the velocity perturbation and the potential perturbation, 
Siu, is not independent of the other variables. For instance, for a single scalar field we have, from Eq. (j4.4()|l . 



S 



lU 



U 



H a^ 



(5.8) 



In a cosmology dominated by the scalar field this forces the scalar field perturbations to become adiabatic, Sju — ^ 0, 
in the fc — > limit for finite \E' But in the presence of other fluids, adiabatic fleld perturbations require the scalar 
field velocity to coincide with the total velocity perturbation, and hence ^^{p-y + P-y)Vij = 0. 

Our formalism is applicable to a variety of cosmological models including interacting scalar fields and fluids. Al- 
though we have only considered minimally coupled scalar fields our formalism can naturally be extended to non- 
minimally coupled fields. Conformally transforming to the Einstein frame introduces interactions between a scalar 
field and the matter, of the form [s^l 

Qf^) oc T„^V''^I , (5.9) 

where is the trace of the matter energy-momentum tensor. On the other hand the perturbed decay of an oscillating 
massive scalar field into light particles may be described by the perturbed energy transfer between pressureless matter 
and radiation |^]. 
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APPENDIX A: GAUGE TRANSFORMATIONS AND SOME USEFUL EQUATIONS 

1. Gauge transformations 

The metric tensor, including linear scalar perturbations about a flat FRW background is given by 

_f-l-2^ aB,, \ 

" ^ aS, a" [(1 - 2^)5,j + 2E^,,] ) ' ^'^'> 
A first order coordinate transformation = + (5a;^, where 5x^ = [5t, 5x^], induces a change in the metric tensor 

where £sxi^ denotes the Lie-derivative with respect to 5x^ . 

The scalar functions </>, t/'i B and E change under the transformation as 

4> = (p-k, (A3) 

^ = ^ + -5t, (A4) 
a 

aB = aB - a^5x + 5t , (A5) 
E = E-Sx. (A6) 
The total four velocity is subject to the constraint 

u^u^^^l, (A7) 

and allowing for linear perturbations we find 

= [-(l + 0),a(« + B),»] , (A8) 



1-0,-1; 
a ' 



(A9) 



A first order coordinate transformation induces a change in the four velocity according to 

Ufj. = Uf, - £sx^Ufi . (AlO) 

In a similar fashion we can define fluid four velocities for the individual fluids. We then find that the velocity potentials 
transform as 

V = V + aSx , Va = Va + aSx , (All) 

and hence the covariant scalar velocity perturbations, defined above in Eqs. (|2.16|) and H2.20|l . as 

Va = a {va + B) , 
V = a{v + B) , 

transform as 

\r = Vc,+St, (A12) 
V = V + 6t. (A13) 

Note, in we used momentum perturbations Sq and Sqa instead of the velocity perturbations V and Va which 
are related by 

Sqa = {Pa+Pa)Va, (A14) 

Sq = {p + P)V. (A15) 
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The total and the individual momenta change as 

6q = 5q+ {p + P)5t, Sqa = Sqa + {pc + Pa)St . (A16) 
The shear scalar, as = a^E — aB, changes as 

(Ts= as - St. (A17) 

The density perturbation changes as 

Sp = Sp-p6t. (A18) 
The perturbed energy transfer four vector of the a-fluid, Q'^^y defined in Eq. (|2.13|l . transforms as 

Q{a) ^ Q(a) - ■^SxfQIa) , (A19) 

and we find that the perturbed energy transfer to the a-fluid changes as 

- 5Q^ - Q^St , (A20) 
whereas the perturbed momentum transfer is gauge-invariant 

7^^f^■ (A21) 

2. Identities 

a. Background 

A useful relation between the Hubble parameter and the background energy density can be derived by combining 
Eqs. (E21), ^M: and to give 

6. Perturbations 

In order to derive the set of evolution equations for Sap and Vap, Eqs. (|3.3HI and (|3.32(l . we made use of the 
following relations between the a-fiuid and total density perturbation and the relative entropy perturbation, 

Pa P iH ^ p 

and the relation between the relative velocity perturbation and the a-fluid and total velocity perturbation. 

In terms of the curvature perturbations on uniform density and comoving hypersurfaces the above expressions give 

= C+^E-'5"7, (A25) 

7 P 

TZa = n-Hy^P^±^Vap. (A26) 
V ^ + ^ 

We found it useful to "symmetrise" some expressions involving differences in the Vap'- 

E - E ^-JT^y^-' = ^ + + ^ (^a - E ^-J^ (^"7 + Vp,) . (A27) 

7 7 7 

A similar expression holds for Safi^ replacing the p^ + P~f and Vap by p-y and Sap- 
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APPENDIX B: EVOLUTION EQUATIONS IN THE UNIFORM CURVATURE GAUGE 

Uniform-density or uniform- field gauges, used to define tlie curvature perturbations Ca and TZi, can become ill- 
defined in some cases even though nothing singular is happening in other gauges. This occurs whenever the uniform 
field hypersurfaces become singular, i.e. if (pi = 0, as can be seen from Eq. 14.25|l |35l l36j |. or when the uniform 
density hypersurfaces become singular, which happens when pa = 0, as can be seen from Eq. H3.1|l . In the latter 
context the equations presented below were used in |3^ . but only in their large scale limit. The singular behaviour 
can be eliminated by simply working with rescaled perturbation variables which physically corresponds to wor king 
with perturbations in the uniform-curvature gauge Ql • The density contrast Ag defined in and denoted A in [3l| 
is related to the density perturbation on flat slices by Ag = A = 6p/p. 

Note that all the variables in this section are defined on uniform curvature hypersurfaces. To avoid confusion we 
denote quantities evaluated in the uniform curvature gauge by a "tilde". The curvature perturbations (, (a, and 
TZa are simply related to the density and velocity perturbations in the uniform-curvature gauge via Eqs. (|3.11(l and 

In the case of multiple fluids the energy evolution equation of the a-fluid in the uniform curvature gauge is given 
from Eq. H2.15|l by 



6po 



3H{l + cl)-f 

Pa 



Spa - SQir, 



5Pw 



The evolution equation for the velocity perturbation of the a-fluid, Va, is, using Eq. H2.17|l 



(1 + 4; 



Pa+p, 



SHci 



Pa 



Pa 



5Pr 



2 

clSpa - o— n 

3 



0. 



(Bl) 



(B2) 



The total velocity and the velocities of the individual species are related by Eq. (|2.19|l . whereas the total density 
perturbation is just the sum of the density perturbations of the individual fluids, Eq. H2.18|l . Summing over the 
evolution equations of the individual fluids, Eq. HB1|I . and using the constraints Eqs. (|2.18|l and (|2.19|) . we find the 

total energy evolution equation to be 



dp + 3H{1 + cl)5p + 3i/^P„ad - ^ {p + P){V + a,) 



0, 



and the total momentum evolution equation from Eq. I)B2I 

' H 



V 



H 



3Hci V 



p + P 



5p 



P 



5P 



nad 



2fc2 
3^ 



rB 



(B3) 



(B4) 



The shear in the uniform curvature gauge, a, is simply the rescaled metric potential '5, defined in Eq. H3.8II . as can 
be seen from Eq. H3.11I) . so to close the system of equations we can cither use the constraint equation 1)3. 7|l together 
with Eq. (|3.11|) . or use the shear evolution equation in the uniform curvature gauge, 



Has 



H 



-V"-87rGn = 0. 



(B5) 
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